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We derive new general expressions for tlie fluctuating electromagnetic field outside a homogeneous 
material surface. The analysis is based on general results from the thermodynamics of irreversible 
processes, and requires no consideration of the material interior, as it only uses knowledge of the 
reflection amplitudes for its surface. Therefore, our results are valid for all homogeneous surfaces, 
including layered systems and metamaterials, at all temperatures. In particular, we obtain new 
formulae for the near-field region, which are important for interpreting the numerous current ex- 
periments probing proximity effects for macroscopic and/or microscopic bodies separated by small 
empty gaps. By use of Onsager's reciprocity relations, we obtain also the general symmetry prop- 
erties that must be satisfied by the reflection matrix of any material. 
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I. INTRODUCTION 

The study of electromagnetic (e.ni.) fluctuations has 
received much attention over the years, because of its fun- 
damental importance in diverse areas of physics. While 
the classical problem of determining the thermal spec- 
trum of e.m. fluctuations in a black body is at the origin 
of quantum theory, the concept of zero-point e.m. quan- 
tum fluctuations leads to new quantum phenomena, like 
the Lamb shift, and the existence of van der Waals forces 
between atoms and/or macroscopic bodies. Electromag- 
netic fluctuation-driven forces become appreciable in the 
submicron range and rapidly increase in the nanometer 
scale. 

It should be observed that while general thermody- 
namical theorems (Kirchhoff's law) give much informa- 
tion on the radiation field far from the body surface, no 
analogously general results are available so far, on the 
problem of characterizing the e.m. field very close to a 
surface. A detailed study of the near-field region is of 
great importance in numerous physical phenomena that 
depend on proximity effects between macroscopic and/or 
microscopic bodies separated by small empty gaps. Im- 
portant examples of such phenomena are the Casimir ef- 
fect [H , radiative heat transfer between closely space bod- 
ies at different temperatures 0, H, 01 , and vacuum friction 

. Non contact interactions mediated by fluctuations of 
the e.m. field are now finding also exciting applications 
to nanotechnology to Bose-Einsten condensates Q, 
and to ultracold atoms [1| . 

The theory of e.m. fluctuations and thermal radia- 
tion was developed long ago by Rytov [2]. In this the- 
ory, the physical origin of e.m. fluctuations resides in 
the microscopic fluctuating currents, which are present 
in the interior of any absorbing medium, according to the 
fluctuation-dissipation theorem. These random currents. 



when introduced into Maxwell's equations, constitute the 
source of fluctuating e.m. flelds, that extend beyond the 
boundaries of the medium, partly in the form of propa- 
gating waves (PW) and partly as evanescent waves (EW) 
localized near the surface. Rytov's method was used by 
Lifshitz [3] in his famous work of van-der Waals attrac- 
tions between macroscopic solid bodies, which consti- 
tutes the basis of all current studies of the Casimir effect 
between real media. The same ideas were also applied 
by Polder and van Hove 3 to investigate the problem of 
radiative heat transfer between closely spaced dielectric 
bodies. While appealing for its conceptual simplicity, 
a drawback of Rytov's approach is that it requires de- 
tailed knowledge of the e.m. fields in the body's interior. 
In particular, the theory requires as inputs the dielec- 
tric tensor ey (aj,x, x') and/or the magnetic permittivity 
liij{u!,x,x') of the material, together with the boundary 
conditions satisfied by the electric and magnetic fields at 
the interface. Therefore, with this approach, the problem 
arises of understanding to what extent the final expres- 
sions for the fluctuating e.m. field outside the body de- 
pend on the assumptions made for the permittivities and 
the chosen boundary conditions, which is a very delicate 
issue in the case, say, of spatially non-local materials. 

The spirit of our approach is radically different, as we 
never consider the e.m. field or the microscopic fluctuat- 
ing currents in the interior of the material body. There- 
fore, we do not make any use of either the constitutive 
equations for the material, or the boundary conditions 
for the e.m. field across the interface. We regard the 
surface as fully characterized in terms of a set of reflec- 
tion amplitudes ~ rxfj^{uj, kj^), depending on the fre- 
quency Lo and the projection k_L of the wave-vector onto 
the surface. These amplitudes are considered by us as 
data, obtained either by experiments or by an indepen- 
dent computation, which is of no concern to us. We shall 
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see that knowledge of the reflection amplitudes is suffi- 
cient to determine, in a completely general way, the corre- 
lators for the e.m. field outside the body. The possibility 
of expressing the correlators of the e.m. field outside a 
material surface in terms of its reflection amplitudes is 
not entirely new. It is clearly suggested by Lifshitz for- 
mula 9], as in its final form it expresses the interaction 
between two nearby parallel dielectric slabs in terms of 
their respective reflection amplitudes. The general valid- 
ity of Lifshitz formula, once expressed in terms of reflec- 
tion amplitudes, has been recognized by several authors, 
and it has been recently demonstrated using the scat- 
tering approach [l^, [lH . The idea that fleld correlators 
can be given in terms of reflection amplitudes was also 
utilized in the study of radiative heat transfer, includ- 
ing evanescent contributions, in Ref. 4]. The distinctive 
feature of our approach is that it relies on general results 
from the thermodynamics of irreversible processes, and 
requires no consideration about the field in the interior 
of the body. The process of reflection is considered as 
a macroscopic irreversible process, entirely described in 
terms of the reflection amplitudes at the surface. There- 
fore, our results are of a very general nature, and they 
apply to arbitrary homogeneous surfaces at all tempera- 
tures, thus providing a generalization of previously avail- 
able results, that includes the case of non-diagonal reflec- 
tion matrices, i.e. mixing TE and TM modes. While for 
PW we recover the well known Kirchhoff's formula, Eq. 
(fTTjl below, we obtain a new general and simple formula 
for the EW contribution, Eq. (fT2|) , which constitutes the 
main result of this paper. By exploiting Onsager's reci- 
procity thermodynamical relations, we obtain also new 
general symmetry properties that must be satisfied by 
the reflection matrix of any material. 



II. THE E.M. FIELD OUTSIDE A 
HOMOGENEOUS SURFACE 



In this Section we derive the expression for the cor- 
relators of the e.m. fleld outside a homogeneous sur- 
face. The procedure we follow is a generalization of the 
method used in Ref. to obtain the Planck radiation 
law. In that paper the authors considered an electric 
dipole in thermal equilibrium with the radiation field of 
a black body. Apart from the intra-dipole binding force, 
the dipole is subjected to two forces of e.m. origin: one 
of them is a systematic damping force, representing the 
reaction to the radiation emitted by the dipole itself, and 
the other is a random force due to the fluctuating electric 
field of the black body (to linear order, magnetic interac- 
tions can be neglected). The damping force is computed 
using Maxwell's Equations, and then one can use the gen- 
eral fluctuation-dissipation theorem to obtain the power 
spectrum of the fluctuating e.m. field, which gives the 
Planck's radiation law. In our case, we are interested 
in determining how the presence of a surface influences 
the radiation field in the vacuum region bounded by the 
surface (for simplicity we limit our analysis to fiat homo- 
geneous surfaces) , and we shall see that the problem can 
be solved in full generality by means of a similar proce- 
dure, but considering, instead of one, two electric dipoles 
placed on a plane parallel to surface. 

To be definite, let us suppose that the body occupies 
the z > half-space, such that its surface is at z = 0. 
By invariance under translations in the x, y plane, the 
fiuctuating electric field in the vacuum region outside the 
body (i.e. in the z < half-space) can be written in 
general as: 



E(«"^Hx^,^) = /°°^ ||^(E(-)(z;c.,kJ+E(^)(z;^,kJ) 



(1) 



where 



E(~)(z;c^,kj.) = [a(°°)(u;,ki)ex + a^^HLu,k^)-{kxe^) 



fai°°)(c^,k^) [r,,e^ + r^.-lk^xe^)! + a(°°) (c^, k^) kp-(kW xe^) + 



E(^)(z;c^,kx) = af^(^,kx)ei + a^^^(c^, kx)-(k^'"^ xe^) 



-.(S)i 



} 6-''=^^ , (2) 
(3) 



In the above Equations, k = k_L -I- /c^z and k'^''^ = propagating towards and away from the surface, = 
k_L — fczZ denote, respectively, the wave-vectors for fields \Jlli'^ jc? — k'j_, and we take the square root such that 
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fc^ > 0, /c" > (with prime and double prime denot- 
ing real and imaginary parts, respectively). We recall 
that real kz correspond to PW waves, while imaginary 
kz describe EW. Moreover, xj_ — xii + yy, e± = zxk^, 
the subscripts s and p denote TE and TM polarizations, 
respectively, and rx^^ ~ r\fi{uj,'kj_) are the reflection am- 
plitudes for an incident field with polarization /i to be 
reflected as a fleld with polarization A. The field E^°°' 
represents an incoming field from z = —oo, that gets re- 
flected by the surface, while we can think of E'^"^-' as the 



fleld radiated by the surface. Since the expansion for the 
magnetic fleld is obtained from Eqs. by integrating 

Maxwell's equations (B — —c J dtVxE), the fluctuat- 
ing e.m. fleld is fully described by the correlators among 

We let a(") , a = oo, 5 

row vector with elements (ai""**, Op"''*). From invariance 
under time translations, and by homogeneity in the x,y 
plane, the correlators must be of the form: 



the amplitudes a^^^"^' (w, kj_) 

a column vector with elements (ai"'' , a^"^ ) , and a^"^^ a 
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i(")(c^,k^),a(^)t(c.',kl)) = (2^)3 C("'^)(a;,k^) Situ ~ lu') S<^'\k^ - k'^) , 

I 



(4) 



with all other correlators vanishing. Therefore, the prob- 
lem reduces to computing the four matrices C^°'^^{uj, k^). 
To do this, we imagine placing in the radiation fleld two 
electric dipoles. We let g^"^-' and A = 1,2 their re- 
spective charges and displacements, and we assume that 
they are placed at arbitrary points x|,^^ and x^-* on a 
plane z — w < in the empty region below the surface. 



For our purposes, it is convenient to suppose that the 
displacements of both dipoles are bound to occur in 
the x,y plane, = Cj^^- For sufflciently small q'^^^ 

(A) 

and , the Equations of motion for the two oscillators 
are linear, and have the following form: 



„ . , -OO 



(^)(i') = g(-')E(f"^)(t,x(L^\H- (5) 



S = l,2 



In this Equation, k'-^^ are the elastic constants for the 
restoring intra-dipolar forces, while the non-local tensor 

t-^(ret) (A.) ( H) 

kernel Ej^ {t — t' , Xj_ — Xj_ , w) represents projection 
onto the x — y plane of the retarded electric fleld pro- 
duced by dipole B at the position of dipole A, which 
produces a systematic frictional force. Finally, the force 
on the r.h.s. is a random term arising from the fluctu- 



ating electric field. The expression of the time Fourier 



transform of the kernel 



,(rot) 



w) 



found by resolving Maxwell Equations in the empty re- 
gion outside the body, a task that is easily accomplished 
by exploiting translational invariance along the surface. 
Upon performing a space-Fourier transform of the fields 
in the z-plane, we obtain: 



J 



E^ (a;;xi^)-xV^',H = ~- , ^^^^ 



27r /■ dki 



+ 



-^(1 + rss e 



— (-1 -f Tppg" 



- ) k^ + rspc - ej_ 



1 e_L 



, (A) (B)s 

,1^^ J, g«k^.(xV-xV) (6) 



r 



Now, Eq. (I5|) is a generalized Langevin Equation, of the 
type commonly used in the theory of irreversible pro- 
cesses [l^l- At thermal equilibrium, it implies a set 
of general relations, collectively known as fluctuation- 
dissipation theorems involving the frictional forces 



on one side, and correlators of random forces or of dy- 
namical variables of the system, on the other. Of par- 
ticular interest to us is the so-called second fluctuation- 
dissipation theorem [l3| . which in our case reduces to 
the following relation between the correlator of the ran- 
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dom force and the time Fourier transform of the damping force: 
I 



Jo 



where F[uj,T) is the quantity 

i^(-,T). -coth(^j ' 

with T the temperature, and ks Bohzman's constant. 
We note that Eq. ([7]) takes account of quantum effects, 
and includes the contribution from zero-point fluctua- 
tions. The expressions for the quantities C^"''^ (w, kx) 
can be now determined by inserting Eqs. (ITJ31) into the 
l.h.s of Eq. (I7|), and requiring that the two members are 
equal for arbitrary values of x^"^^ , x^^^^ and w. After some 
tedious but straightforward computations one obtains for 
C("'3)(w,kj_) the following expressions: 

)=n^,r)^Rcfl) , (9) 



Eq. ^ coincides with the emission spectrum of a black 
surface, and Eqs. pO)) show that this radiation is uncor- 
related with that emitted by the surface. As for C^^^\ 
when kz is real, i.e. for PW, we find: 

^ ^j^(^^y)(l_^^t) ^ (PW) (11) 
while for kz imaginary, i.e. for EW, we obtain: 

C^'''^^-^^,F{u:,T){R-R^), (EW) (12) 

C \Kz\ 

where R is the reflection matrix i?^^ — r\f^. Eqs. (|lip 
and (|12p describe the emission from the surface. It should 
be observed that the contribution from PW, Eg. pTj) is 
equivalent to Kirchhoff 's law, and it is obviously possible 
to derive it by simply imposing detailed energy balance 
between the radiation flux from the black surface at in- 
finity and the flux radiated by the surface. The same 
type of argument, however, is useless to determine the 
EW contribution, Eq. p^ . because the associated aver- 
age flux of energy is zero. Therefore, Eq. (|12p should be 
regarded as a non-trivial extension of Kirchhoff 's formula 
to the near-field region, and Eqs. PT|) and together 
provide the complete description of the e.m. field outside 
the body. It should be noted that EW fluctuations are 
zero when the reflection matrix R is hermitean. 

We remark an interesting consequence of the symme- 
try properties of the correlation functions, in Eq. ([7]). In 



the absence of applied external magnetic fields, the cor- 
relators on the l.h.s. of Eq. ([7]) are invariant under time 

reversal, and therefore the quantities Eij on the r.h.s. 
satisfy the following symmetry property: 

Eij (uj;x\' -x\\w} ^Eji (w;x^ ' - x^ % w) . 

(13) 

These relations constitute an example of Onsager's reci- 
procity relations. Upon using Eq. ([5]), and recalling that 
Eq. must be satisfied for all x'^'^^x^^^ and w, one 
finds that Eq. ([13]) are equivalent to the following sym- 
metry conditions for the reflection amplitudes: 

rssiu!,kj_) = rssiui,-k±) , (14) 
rpp{uj,k±_) = rpp{uj,-k^), (15) 
rsp{uj,k±) = -rps{uj,-kj_) . (16) 

The above relations can be used a criterion to reject spe- 
cific models for the reflection amplitudes. An interesting 
example of this sort is provided by the Drude-Born model 
for chiral materials, which is based on the following con- 
stitutive Equations (in the frequency domain): 

D = eE-/VxE, B = H. (17) 

It can be shown that the reflection amplitudes implied 
by this model do not satisfy Eqs. and therefore this 

model must be rejected. However, the model for chiral 
materials proposed by Fedorov [ij] 

D = e(E + /3VxE) , B = ^(H + /3VxH) (18) 

satisfies the above conditions, and is therefore viable 
from the point of view of thermodynamics. Moreover, 
it is worth noting that chiral Fedorov's materials, even if 
transparent, have a non hermitean reflection matrix R. 
Proximity effects related to EW near the surface of chiral 
media are therefore expected. 

III. CONCLUSIONS AND DISCUSSION 

In conclusion, using general arguments from the ther- 
modynamics of irreversible processes, we have obtained 
general formulae for the fluctuating e.m. fields outside a 
homogeneous surface. Our results are valid for all ma- 
terials, at all temperatures, and provide a generalization 
of Kirchoff's law to the near field region. In particular. 
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they apply to non-diagonal reflection matrices, that mix 
TE and TM modes. Moreover, by use of Onsager's reci- 
procity relations, we have also obtained the general sym- 
metry properties that reflection amplitudes must satisfy, 
in order to be consistent with the requirements of ther- 



modynamics. The same methods can be also applied to 
cavities bounded by homogeneous surfaces, and to sys- 
tems out of thermal equilibrium, but we postpone to a 
separate paper the discussion of these problems. 
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